function s = qam(data,symbolRate,sampleFreq,M)
                
DataConstel = qammod(data, M, 'Gray')
NperBytes=round(sampleFreq/symbolRate);
s=repelem(DataConstel,NperBytes);

% N = length(DataConstel);
% K = sampleFreq/symbolRate;
% s=zeros(1,N*K);
% for i=1:K
%     s(i:K:end)=DataConstel;
% end